The Hubbard model as an approximation to the entanglement in nanostructures 
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We investigate how well the one-dimensional Hubbard model describes the entanglement of particles trapped 
in a string of quantum wells. We calculate the average single-site entanglement for two particles interacting via 
a contact interaction and consider the effect of varying the interaction strength and the interwell distance. We 
compare the results with the ones obtained within the one-dimensional Hubbard model with on-site interaction. 
We suggest an upper bound for the average single-site entanglement for two electrons in M wells and discuss 
analytical limits for very large repulsive and attractive interactions. We investigate how the interplay between 
interaction and potential shape in the quantum well system dictates the position and size of the entanglement 
maxima and the agreement with the theoretical limits. Finally we calculate the spatial entanglement for the 
quantum well system and compare it to its average single-site entanglement. 

PACS numbers: 03.65.Ud, 71.10.Fd, 73.21.La, 73.21.Fg 



I. INTRODUCTION 

Entanglement is considered one of the main resources in 
quantum information and a reason why quantum computers 
may be used for computing feats that could not be achieved 
with traditional processors In this context quantum dots 
are thought of as a viable possibility in the quest to construct 
scalable quantum processors |2-9]. With this in mind, finding 
accurate ways to calculate the entanglement between electrons 
in quantum dots becomes important for quantum information 
processing. However modeling these many-body quantum 
systems often necessitates the employment of approximations. 
For example, one-dimensional wells may be used in the study 
of spherically symmetric quantum dots and aid understanding 
of trends in more general quantum dot systems 1 10]. 

The Hubbard model [II] allows interacting many -body sys- 
tems to be simulated by mapping them onto a lattice model 
with (usually) on-site interactions only. Despite its relative 
simplicity it still captures a significant amount of physics: 
for example, in strongly coiTelated fermionic systems it has 
been used to model particles trapped in an optical lattice 
highTc superconductivity [13], and the metal-insulator transi- 
tion [14]. The one-dimensional homogeneous Hubbard model 
(HM) also benefits from the existence of an exact solution in 
the thermodynamic limit ifTsll . Recently the use of the Hub- 
bard model as an approximation to the exchange coupling 
in quantum-dot nanostructures has been investigated [161 . 
The entanglement of the one-dimensional Hubbard model has 
been investigated in Refs. iflTl - fTgll . A local-density approx- 
imation (LDA) to the entanglement has been proposed in 
Ref. ||2^ and applied to inhomogeneous systems. 

In this paper we compare the Hubbard model predictions 
to results from a system of two interacting fermions trapped 
within a chain of square well potentials. Here each well cor- 
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responds to one of the Hubbard model sites. We consider both 
repulsive (electron-electron) and attractive interactions, and 
calculate the coiTesponding entanglement when the strength 
of the interaction, the chain length, and the interwell dis- 
tance are changed. We compare these results with the average 
single-site entanglement calculated from the Hubbard model. 



In doing this we infer information on the accuracy of using 
results from the Hubbard model to approximate the average 
single-site entanglement of the quantum well system. If the 
entanglement of electrons in quantum wells can be described 
using the Hubbard model, then, by using the powerful LDA 
formalism developed in ll20ll . we could in principle calculate 
the entanglement present in quantum well systems with a large 
number of interacting electrons. This would be a non-trivial 
result as a direct calculation of the entanglement for a system 
with a large number of interacting electrons becomes compu- 
tationally prohibitive as the number of particles increases. 



In Sec.|II]we introduce the one-dimensional quantum well 
system and discuss how we numerically calculate the single- 
site entanglement. In Sec.|III]we compare the results from the 
Hubbard model with the ones from the electron system for dif- 
ferent interwell distances, chain length, and Coulomb interac- 
tion strength. We also investigate how the interaction strength 
affects the electronic density, and its effect on the matching 
between the Hubbard and quantum well system. In SeclIVI 
we propose an upper bound for the average single-site en- 
tanglement of the quantum well system and discuss the large 
Coulomb interaction limit. We investigate how close our nu- 
merical results come to these analytical expressions. In Sec.lVl 
we compare attractive and repulsive particle-particle interac- 
tion, as well as discussing the large interparticle attraction 
limit. Sec. |VT]is devoted to the comparison between the av- 
erage single-site entanglement with the spatial entanglement 
for the quantum well system and finally Sec lVIII contains our 
conclusions. 
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II. THE QUANTUM WELL TWO-ELECTRON SYSTEM 

The Hamiltonian, in atomic units, for the one-dimensional 
system of two electrons confined within an array of quantum 
wells (QWs) is given by 



H 



i=l,2 



2dxi 



The potential v{xi) models a string of regularly spaced, iden- 
tical square wells, symmetric around the origin, and defined 
by the quantities: M the number of wells, w the width of each 
well, d the barrier width between two consecutive wells, and 
Wo the depth of each well. We set f{x) — 5{x) to model a 
contact Coulomb repulsion and define Cu as the interaction 
strength. Cu will allow us to compare the system with the 
Hubbard model. 

We solve the time independent Schrodinger equation corre- 
sponding to Eq. ^ by using 'exact' diagonalization; the elec- 
tronic ground state is a singlet thereby satisfying the stipula- 
tion of zero magnetization. We calculate the spatial part of the 
many-body wavefunction using the first N eigenfunctions of 
the potential v{x) as single-particle basis functions. We em- 
ploy these to produce a basis of symmetric two-particle wave- 
functions which means we only need to consider N{N + 1)/2 
functions. The form of the Hamiltonian is conducive to this 
method as by varying the interaction strength independently 
we do not need to recalculate the basis, or any integral in- 
volved in the diagonalization, as could be the case if we varied 
the well geometry directly. In this respect we note that a sys- 
tem with Cu ~ K, depth vq, well width w and barrier width 
d is equivalent to a system with Cu ~ 1, depth vq/K'^, well 
width Kw, and barrier width Kd. 



A. Average single-site entanglement 



with Pi (7) the probability of double (7 ='ti), single (7 =t or 
4,), or zero (7 = 0) electronic occupation at site i lfl9ll . 

We calculate the ground-state wave-function, for an even 
number M of wells, and from that obtain the occupation prob- 
abilities. We calculate the probability that both electrons are 
in the same left-most (Afth) well as 



Pm{U)^I / \^lj{xi,X2)\'^dxidx2, 



— Tc ^ — r, 



(5) 



where Vc is the (numerical) integration cut-off point and h = 
— (M/2 — l){w + d) the mid point between the left-most well 
and the next well. 

The probabihty that only one spin up (or spin down by sym- 
metry) electron is in this well is 



PM{t)^PM{i)^ / / \^ixuX2)rdxidX2. (6) 

Pm{0) may then be deduced, as the probabilities sum to 1. 
For the other wells the occupation probabilities are 



Pj+M/2{tl^) — 



b pb 



and 



(t) 



b fa 



a J b 



\ijj{xi,X2)\'^dxidx2, (7) 



\ip{xi,X2)\^dxidx2 (8) 



\ip{xi,X2)\'^dxidx2, (9) 



with a = —j{'w + d) and h = —j{w + d)+d + w, 
1 < j < {M/2 — 1). As we only consider an even number 
of wells distributed symmetrically about the origin, the prob- 
ability values for wells M/2 to 1 are known by symmetry. 



We wish to calculate the average single-site (or local) entan- 
glement of the system ground state. This type of entanglement 
is relevant for systems of indistinguishable fermions |21]. To 
this aim we divide our QW system in contiguous 'sites', each 
site centered around a single well. 

The entanglement entropy S of the system is given by 

1 

i 

with 

St = ~Trpred,i l0g2 Pred,i (3) 

the /-site von Neumann entropy of the reduced density ma- 
trix Pred,i- The von Neumann entropy is considered as one 
of the definitive measures of entanglement for a pure bipartite 
system. By dividing the system into sites and moving to a site- 
occupation basis the reduced density matrix pred.i becomes a 
4x4 diagonal matrix |21, 22], 

Pred,^ = diag [P,(n), P.(t), Pt{i) . ^^.(0)] (4) 



III. COMPARISON WITH THE HUBBARD MODEL 
The Hubbard model ifTTll is described by the Hamiltonian 

+^7^n,,tf^^4' (10) 

i 

where i runs over the M sites and a —^^ J,- Here t is the hop- 
ping parameter and U is the interaction strength. c| ^ (ci^cr) 
creates (destroys) a particle of spin a at site i while n^ cr = 
c| ^Ci^a- is the particle number operator We solve Eq. (fTOl i 
by exact diagonalization in the single-site occupation basis 
{ Iti) I It) I \i) 1 10)}. We apply open boundary conditions and 
consider an average particle density of ri = + = 2/ M, 
with Ua- the average density of the cr-spin component. Again 
we calculate the average single-site entanglement Eq. (|2]i [ l?]- 
I19II . Usually the hopping parameter t is used to rescale IJ, 
giving the dimensionless interaction strength U — U/t. 
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To compare results, we need to calculate the equivalent of t 
and U for the QW system discussed in Sec. [Ill In the Hubbard 
model the hopping parameter is the expectation value of the 
single-particle operators in the Hamiltonian with respect to 
the single-particle wave functions localized at adjacent sites. 
When the hopping parameter t is independent of the sites it 
may be written as 

t = {Mr)\ (^-^V2 + V{r)^ |0.+i(r-)) , (11) 

where (pi is the wavefunction at any site i and V{r) is the 
single-particle confining potential. 

Following this definition, we estimate the hopping parame- 
ter for our quantum well model as 

where (j)L(R) has the shape of the single-particle ground state 
of the finite single square well potential (p^j, but centered in 
the left ((pL) or right ((p^) well. Here v{x) is the potential de- 
fined in Eq. with the zero of energy chosen such that v{x) 
has zero as its lowest value thereby ensuring that the potential 
contribution is always positive and an increase in the depth of 
the well causes the hopping parameter to decrease. The phase 
of (f>L and (pji is chosen to make real and positive. 

The on-site interaction in the Hubbard model lil ill is defined 

by 

U^l {Mr2)\ {Mri)\ l'^^(^i)) . (13) 

The corresponding parameter in our ID model with a delta 
function interaction is then 

U^^^ j <Pt{x)dx. (14) 
We may now estimate U for our model as 



For systems of wells characterized by the parameters w ~ d ~ 
2 qq, where oq is the (effective) Bohr radius, and vq — 10 (ef- 
fective) Hartree, we find U — 123AACu', when the interwell 
distance is reduced to d = 0.2 oq we obtain U = 3.14Cu, and 
U = 1.27C(7 for the limiting case d = 0. 
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FIG. 1: Upper Panel: electron density for f/ = (Cu ~ 0) (dashed 
line) and f/ = 40 (Cu = 0.00324) (solid line) for a 4 well potential 
and w = 2 ao, d — 2 ao and vq = 10 Hartree. Lower Panel: as for 
the upper panel but for a system of 8 wells. 



electron repulsion, the difference between the electron den- 
sity in the inner and outer wells becomes much smaller (solid 
line). For eight wells and U ~ the inner wells are again 
preferentially occupied while there is very little density in the 
outermost wells (Fig.[Tl lower panel, dashed line). When the 
interaction is 'switched on' to U — 40 the two central wells 
have a lower peak density than the nearby wells to compen- 
sate for the Coulomb repulsion, but the outermost wells still 
display, by comparison, a much lower density. Similar behav- 
iors of the interacting and non-interacting densities are found 
when the distance between wells is reduced by an order of 
magnitude (d — 0.2 ao). However in this case the density is 
considerably different from zero in the barrier region (Fig.©. 

These results seemingly show that, apart from accidental 
compensation, the electron density in the different wells can 
not be made equal by applying an unmodulated Coulomb in- 
teraction. 



B. Comparison of entanglement results 



A. Effect of many-body interactions on the electronic density 

Next we explore how the electron density alters to maxi- 
mize its exposure to the attractive confining potential whilst 
attempting to minimize the interaction between the electrons. 
For two wells the density profile is clearly symmetric so we 
do not discuss it further A non-interacting four well system 
(upper panel of Fig. [T] dashed line) displays a density clearly 
higher in the inner wells. However for U = 40, due to the 



In Fig. Owe compare the average single-site entanglement 
for the Hubbard model with that of the QW electron systems 
characterized hy d = 2, d — 0.2, and the limiting case d = 0. 
The latter corresponds to the arbitrary partition of a single well 
of width Mw into M equal regions. 

For d = 2 ao and two wells (upper panel) we see that the 
entanglement decreases similarly to the Hubbard model as U 
increases but the entanglement in the QW system is slightly 
higher When we consider four wells (middle panel), in both 
cases the entanglement increases up to a maximum and then 
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FIG. 2: Upper panel: electron density for (7 = {Cu ~ 0) (dashed 
line) and U = AO (Cu = 12.75) (solid line) for a 4 well potential 
with w — 2 ao, d — 0.2 ao and vq = 10 Hartree. Lower Panel: as 
for the upper panel but for a system of 8 wells. 



decreases, the maximum occurring at slightly different values 
of U. For eight wells (lower panel) the entanglement in the 
two systems is almost indistinguishable. 

For d = 0.2 ao, two wells, and U < 8 the Hubbard model 
is fairly accurate in reproducing the average single-site en- 
tanglement, while, for stronger interactions, results from the 
Hubbard model reproduce only the qualitative trend (Figl3j 
upper panel). For U > 8 the entanglement values are interme- 
diate between the Hubbard model and the limiting case d — 0: 
in this respect we note that for d — 0.2 uq even when the re- 
pulsive interaction is as high as U — 40 (Cu — 12.75) the 
electron density in the QW system does not become localized 
within the wells (see Fig.|2|i. 

When four wells are considered, the Hubbard model repro- 
duces the entanglement trend qualitatively and is less accurate 
when the interaction is low (Fig. [3] middle panel). The max- 
imum entanglement is lower for d = 0.2 ao and in general 
the entanglement trend is intermediate between the Hubbard 
model and the limiting case d = 0. For four wells the differ- 
ence between the maximum values of the entanglement can 
not be removed by rescaling U (see discussion in next sec- 
tion). 

For eight wells the Hubbard model reproduces the qualita- 
tive behavior but the entanglement is lower at all values of U 
and intermediate with respect to the results for d = 0. 

It should be noted however that, even for d ~ 0.2 ao, the 
percentage error for the entanglement as estimated using the 
Hubbard model will be relatively small for the eight and four 
well systems (~ 1%), while more substantial for the two-well 
case and [/ > 8 (- 20%). 
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FIG. 3: Average single-site entanglement for the Hubbard model 
and the QW electron system with U — Uw/tw, w = 2 ao, d = 2 
ao, d — 0.2 ao or d = 0, and vo — 10 Hartree. Upper panel: 2 sites 
with n = 1 ag^ and 2 wells. Center panel: 4 sites with n = 0.5 a^^ 
and 4 wells. Lower panel: 8 sites with n — 0.25 a'^,^ and 8 wells. 



Our results show that the average single-site entanglement 
of the Hubbard model is a very good match for the entangle- 
ment of a QW electron system when wells are far enough apart 
to prevent significant electron density in the interwell barrier 
region (see Fig. [T]); it is less good, although it gives the gen- 
eral trend, when the wells become closer, as the density profile 
displays less well-defined 'sites' (see Fig. |2]) . Surprisingly 
though, when considering a large number of wells, the Hub- 
bard model reproduces the entanglement within few percent at 
all interaction strengths, even when compared to the limiting 
case scenario d = (FiglSj lower panel). This results sug- 
gests that here the Hubbard model sites could be interpreted 
as a fine enough mesh discretization of the continuous spatial 
variable. 



C. Rescaling (7m /tuj 

We now investigate whether, for d — 2 uq, the small dis- 
crepancy between the Hubbard model and the QW system re- 
sults for the entanglement may be removed by choosing an 
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'ad hoc' value of 11^/ (twCu). 

We find that with U^/ {Cjjtw) = 11500 the entanglement 
for the QW system is almost identical to the results from the 
Hubbard model for all the systems considered (Fig.|4|i. 




FIG. 4: Average single-site entanglement S for the Hubbard model 
and the QW electron system with U = llSOOCiy, w = 2 ao, d = 2 
ao, and vq = 10 Hartree. Upper Panel: 2 sites with n = 1 ag^ and 
2 wells. Center Panel: 4 sites with n = 0.5 a^^ and 4 wells. Lower 
panel: 8 sites with n — 0.25 a^^ and 8 wells. 



This suggests that although the calculated Uw/tw gives a 
good estimate for the parameter U to be used in the Hub- 
bard model, we may improve the entanglement accuracy by 
fitting Uw/iCjjtw) and thereby compensate for some of the 
differences between the models. The extent of the scaling 
confirms that — at least for parameters for which there is no 
significant electron density in the barrier regions — the use of 
a single square well wave-function is a good approximation 
in the calculation of Uw/t^ as the result is very close to the 
scaled value. 



IV. UPPER BOUND FOR THE ENTANGLEMENT AND 
STRONG COULOMB INTERACTION LIMIT 

Let us consider the case of zero magnetization, i.e. P(t) = 
P{i) and 2A/' particles where TV is an integer. With M wells, 
we then have a constraint from conserving the particle number 



M 

E 



(16) 



and constraints from the requirement that occupation proba- 
bilities for any well/site must sum to one 



A = ^'.(t;) + 2P.(t) + p.(o) - 1 = 0. 



(17) 



We use Lagrange multipliers to maximize S subject to these 
constraints, i.e. 



d 



S-X(j) 



M \ 



(18) 



with 7 =t, 'fl, 0. Eliminating A and fii from the resulting 
equations give 



p,in)p^io) = (Em' 



(19) 



Eq. (fT9] l relates occupation probabilities within each site, so 
we can find a local maximum of the entanglement where 
all the wells/sites are equivalent, i.e. Pi{j) — P{j)- Irn- 
posing this condition on Eqs. ( fTSI l. ([TtT i and ( fT9] l gives 
P(t) = M/M - {M/Mf and P(n) = {M/lVlf. We note 
that it is only the ratio Af/M that matters for the probabilities 
and hence the entanglement. For 2 particles, the largest entan- 
glement occurs for two well/sites (Fig. O; this suggests that 
JV/M =1/2 (half-filling) may be the condition to obtain the 
largest maximum for the entanglement. 

We now continue with Af = 1 and note that 
P(1-) = 1/Af - l/A/2 and P{U) = ^/M^ correspond to the 
condition of no preferred well and two uncorrelated particles 
of opposite spin. 

Reduced density matrices with eigenstates equal to 1/g, 
g the number of degrees of freedom, would correspond to 
maximal entanglement. However, under the stipulation of 
preserving the particle number together with the request that 
J2j Piil) = 1' '^his state cannot be achieved except for 
M = 2. We could think of moving closer to this by attempt- 
ing to achieve reduced density matrices with more homogene- 
ity within the eigenvalues, i.e., achieving ^^(7) « Pi (7'), at 
least within certain wells. We implement this by relaxing the 
condition that the wells are equivalent and moving part of the 
particle density from one site to another. We reduce Pi (f) and 
Pi (4,) by q while increasing Pj (^4) at site j ^ i hy q, with 
the empty occupation probabilities adjusted accordingly. Set- 
ting dS/dq = gives q = suggesting that the maximum 
entanglement occurs when all wells are equivalent. 

Under this condition, the average single-site entanglement 
is given by Eq. Q, and simplifies to 



5:^(M) = 21og2(Af) 



^ - 1 ) log2 (Af - 1) . (20) 
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M 


qth 

'-'max 


5'max 


2 


2 


2 


4 


1.623 


1.550 


6 


1.300 


1.234 


8 


1.087 


1.033 



M 


"5(7— >oo 


Su=io 


5(7=450 


2 


1 


1.030 


1.000 


4 


1.5 


1.503 


1.500 


6 


1.252 


1.226 


1.226 


8 


1.061 


1.032 


1.032 



TABLE I: Table showing the maximum theoretical average single- 
site entanglement Eq. i20\ . and the maximum entanglement as calcu- 
lated for the QW electron system for d = 2 and different numbers 
of wells. 



TABLE IIL Table showing the limiting value for the average single- 
site entanglement entropy Eq. illi . and the results from the QW elec- 
tron system with U — and U — 450 for d = 2 uq and different 
numbers of wells. 



This maximum average single-site entanglement decreases as 
the number of wells increase (see Table and, for two parti- 
cles, 5*3^ M-^co Q ^j^^ Hubbard model picture, this would 
correspond to the limit of the number of sites going to infin- 
ity and the average particle density going to zero. This limit 
would in fact be expected to have no entanglement as it is es- 
sentially a product state of empty occupations. 

For d — 2 Go, 5'*3x is reached for U = and two wells, 
similarly to the Hubbard model; however, for M > 2, some 
interaction is required to balance the propensity of the non- 
interacting wave-function to favor inner wells. Turning on the 
repulsion between electrons will tend to reduce the discrep- 
ancies between the electron density peaks in different wells 
(see e.g. the upper panel of Fig. [T]); however this will also 
tend to decrease the double occupation probability and in par- 
ticular the already too low value of Pi(M)(ti) the outer 
wells. Therefore, due to the open boundary conditions here 
considered, it may not be possible for the system to reach the 
theoretical maximum for the entanglement by simply varying 
U, as a perfect balance between occupation probabilities in 
different wells may not be achieved without, for example, a 
spatial modulation of the particle-particle interaction. 



d (ao) 


A(t) 




pth 


0.2 


0.157 


0.277 


0.1875 


2 


0.164 


0.265 


0.1875 


d (ao) 


Piin) 


P2in) 




0.2 


5.62 X 10"^ 


0.0595 


0.0625 


2 


8.58 X 10"=* 


0.0626 


0.0625 


d (ao) 


A(0) 


P2(0) 


PntiO) 


0.2 


0.681 


0.386 


0.5625 


2 


0.663 


0.408 


0.5625 



TABLE IL Occupation probabilities for M = 4, interwell distances 
d — 0.2 ao and d = 2 ao, and interaction value Usma^ correspond- 
ing to the entanglement maximum. Theoretical values as used in 
Eq. ^ 

In Table HI] we compare the occupation probabilities 
P\wdx{l) corresponding to the maximum theoretical entangle- 
ment 5*ax to the occupation probabilities calculated for the 
maximum value of the entanglement for the M = 4 system 
and interwell distances d — 2 and d — 0.2. We note that the 
largest discrepancy with the theoretical values is observed for 



the double occupation probability, with Pi(t4.) < -P*ax(ti)- 
This is greatly responsible for the fact that S'^ax < -Smax for 
both interwell distances. The d — 0.2 system also presents 
the largest discrepancies \Pi{'j) — Pmax{l)\ for ^^y 7 and i, as 
well as a large inhomogeneity between well occupation prob- 
abilities I Pi (7) — P2 (7) I for any 7. These account for the fact 
that this system presents a lower maximum for the entangle- 
ment in respect to the system with d = 2 gq. 

We may calculate the theoretical limit for the entanglement 
of two electrons and M wells when U 00 and all wells are 
equally favorable. Using a similar procedure to section UVj but 
with P{U) = and P(t) = P(4.) = 1/M, we obtain 

^*^oo = -|log2(M) + 

We see in Table [III] that Eq. ( I2TI ) describes the large U limit 
of the QW system fairly well, with a percentage error of at 
most 3%. For M = 6 and M = 8 the entanglement of the 
QW system saturates at L/ w 40 and remains slightly below 
the theoretical limiting value as in this case the assumption of 
equivalent wells does not hold even for very strong interac- 
tions (Fig.[T] lower panel). 

In Table |IV] we explore the differences between the theo- 
retical limiting results and the M = 4 system. We consider 
d = 0, d ~ 0.2 ao, d — 2 ao, and the results from the Hub- 
bard model. Table IIVI shows that the occupation probabilities 
for d = 2 ao are almost identical to the Hubbard model and 
extremely close to the theoretical limiting values. For d — 0.2 
instead, no matter how strong the Coulomb repulsion between 
particles is made (U ~ 450 in the table), the very narrow in- 
terwell barriers fail to counteract the effect of the boundary 
conditions, which favor occupation in the central wells. In 
general the inhomogeneity between well occupation probabil- 
ities I Pi (7) — P2(7)| increases for decreasing d, underlining 
the fact that the definition of 'sites' become more arbitrary. 
However the substantially larger double occupancy probabil- 
ity encountered for d < 2 increases the available degrees of 
freedom and hence the entanglement. This confirms the en- 
tanglement trend observed in Fig. |3] center panel. 
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d (ao) 


A(t) 




P*'^(t) 





0.196 


0.284 


0.25 


0.2 


0.228 


0.268 


0.25 


2 


0.249 


0.251 


0.25 


Hubbard 


0.249 


0.251 


0.25 


d (ao) 


A(ti) 


P2{n) 







2.04 X 10"" 


1.81 X 10~^ 





0.2 


2.12 X 10"'' 


3.46 X lO"'* 





2 


6.80 X 10"'^ 


5.87 X 10"^ 





Hubbard 


4.94 X lO"'^ 


4.43 X lO"'* 





d (ao) 


A(0) 


P2(0) 


P"'(0) 





0.608 


0.414 


0.5 


0.2 


0.544 


0.461 


0.5 


2 


0.501 


0.499 


0.5 


Hubbard 


0.501 


0.499 


0.5 



TABLE IV: Occupation probabilities for M = 4, interwell distances 
d = 0, d = 0.2 ao and d = 2 ao, and U — 450. The correspond- 
ing values for the Hubbard model are reported as well. Theoretical 
values as used for Eg. 12 II 



V. ATTRACTIVE VERSUS REPULSIVE 
PARTICLE-PARTICLE INTERACTION 

We wish to discuss how the entanglement pattern is modi- 
fied when we compare attractive {Cjj, U < 0) with repulsive 
particle-particle interaction. In the following we will consider 
the QW system with d = 2 oq and the Hubbard model. In 
Fig. |5] we show the change in the average single-site entan- 
glement with U for different numbers of wells. From our 
calculations Snax always occurs for U > and corresponds 
to [/ = for two, U = 2.1 for four, while U = 4.8 for 
eight wells. As expected from Eq. (|20] |. the maximum aver- 
age single-site entanglement 5max decreases with increasing 
number of wells and our conjectured theoretical maximum 
entanglement S^^^ is indeed an upper bound, to which the 
actual system comes reasonably close (Table For M > 2, 
due to the non-periodic nature of the system, an unmodulated 
interaction strength drives the system towards having equiv- 
alent wells only in the very large \U\ limit. However, as 
particle-particle interaction would naturally introduce corre- 
lations, any spatial modulation of U should probably be non 
trivial in order to mimic the uncorrected electrons' occupa- 
tion probabilities corresponding to the maximum theoretical 
entanglement Eq. (|20] |. 

For U < the entanglement decreases monotonically for 
increasing \U\. This is due to a disproportionate increase of 
the double occupation probabilities, which are favored by the 
attractive interparticle interaction. This limits the access to 
other degrees of freedom which might contribute to the entan- 
glement, and consequently the entanglement is reduced. 

Our calculations show that the Hubbard model reproduces 
well the average single-site entanglement of a QW system 
with relatively wide interwell barriers. The comparison for 



M 


'-'!7<<0 


c'h,2 
^U<<0 


S'!7=-40 


2 


1 


1 


1.030 


4 


0.5 


0.811 


0.706 


6 


0.333 


0.65 


0.585 


8 


0.25 


0.544 


0.497 



TABLE V: Table showing the theoretical limits for the average 
single-site entanglement and U <^ for different numbers of wells. 
Results from the QW system with d = 2 ao and U = —40 are 
presented as well. 



M = 4 is shown in Fig. |5] 





/ '\ QW system 4 wells' 




Hubbard 4 sites 




/ \ QW system 2 wells 




\QW system 8 wells 





'■-40 -20 20 40 

U 



FIG. 5: Average single-site entanglement of the 4-sites Hubbard 
model and of the QW system with 2, 4, and 8 wells vs U , U = 
Uw/tw, d = 2 ao, w — 2 ao, and no = 10 Hartree. 



A. Large inter-particle attraction limit 

For U << the two center wells could have equal proba- 
bilities of double occupation and emptiness whilst all the other 
wells would be empty. This would lead to an average single- 
site entanglement of 

= 2/M. (22) 

We see in TablelVlthat the QW system does not get very close 
to this limit except for M — 2. The form of the confining po- 
tential is such that all the wells will always contain some den- 
sity for the finite interaction strengths considered {U > —40). 
At these interaction strengths the system is better described 
by assuming that all wells are equivalent but that there is no 
single occupation. This gives 

^*'<<o = ^ log. (M)- (l - ^) log, (l - ^) . (23) 

Fig. |5] shows that the entanglement remains intermediate be- 
tween Sf}'^^f^ and 5'*<<0' '^^^ relatively Umited effect 
of the short-range interaction considered. 
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VI. SPATIAL VERSUS AVERAGE SINGLE-SITE 
ENTANGLEMENT 

In this section we consider a different type of entangle- 
ment contained within the QW system, the spatial entangle- 
ment between the two trapped particles. This represents the 
particle-particle entanglement spanning from the many-body 
wave-function spatial degrees of freedom. Once more we cal- 
culate the entanglement using the von Neumann entropy of 
the reduced density matrix, Sgp ~ — 2^?'Pred,sp log2 Pred,,sp, with 
Pred.sp calculated from the spatial degrees of freedom as l(23ll 

Pred,sp(a;i,2:2) = j [xi.Xy,)'^ {x2,X'i)dX'i- (24) 

This expression is diagonalized with respect to the basis set 
employed. Also in this case we allow for attractive as well as 
repulsive interaction between the particles. 

Notice that in the present case the spatial entanglement is 
zero when there is no interaction as the wave-function fac- 
torizes into spatial and spin components and the implicit cor- 
relations arising from the Pauli exclusion principle — and the 
related entanglement — are accounted for within the spin de- 
grees of freedom. 

For two wells (Fig. |6] upper panel) we see that the spatial 
entanglement is a mirror image of the site entanglement when 
reflected along the line 5, Sgp = 1. For larger numbers of 
wells the relationship is more complicated. In most regions 
when the spatial entanglement increases the site entanglement 
decreases and vice versa; however the spatial entanglement 
does not have a minimum exactly where the average single- 
site entanglement has a maximum. This is because the spatial 
entanglement's minimum always occurs at = when there 
is no correlation between the particles' positions. 

An intuitive explanation for the almost opposite behavior of 
these two types of entanglement is that for [/ > increasing 
U increases the repulsion and the correlation between parti- 
cles. Hence one electron's position reveals more about the 
other electron while the number of spatial degrees of free- 
dom is not dramatically limited, so the spatial entanglement 
increases. However the probability of double occupation is 
reduced by a large positive interaction so less is learned by 
measurement with respect to wells/sites even though the elec- 
tron affects the other's position more. Therefore the site en- 
tanglement decreases once it has reached its maximum but, 
for M > 2, much less strongly than the spatial entanglement's 
increase. 

For U < Q the reduction in the probability of single occu- 
pation causes the average single-site entanglement to decrease 
markedly when \U\ increases. The increase in spatial entan- 
glement with increasing \U\ here comes from the system ap- 
proaching the situation where measurement of one electron's 
position reveals the other electron to be in the same region. 
This results in large entanglement and we find that the spa- 
tial entanglement for U < Q increases as the number of wells 
where the electron can be found increases. 










Swells 




Average single site 




Spatial 









"^^20 ' -10 ' ' 10 ' 20 



U 

FIG. 6: Average single-site entanglement S and spatial entangle- 
ment Ssp for the QW electron system with U — Uw/tw, d — 2 
ao, w = 2 ao, vo = 10 Hartree and 2 wells (upper panel), 4 wells 
(center panel), and 8 wells (lower panel) 



VII. CONCLUSION 

In this paper we examined the average single-site and spa- 
tial entanglement of two particles confined in a string of quan- 
tum wells and interacting via a contact interaction. The re- 
sults for average single-site entanglement were compared to 
those of the one-dimensional Hubbard model with on-site in- 
teraction, to investigate when this model is a good approxi- 
mation to the two-particle system. For repulsive (Coulomb) 
interaction, we found that the trend of the entanglement was 
reproduced, with a generally good quantitative agreement, 
when comparing with a Hubbard model characterized hy U = 
Uw/tw where C/^u and were calculated from the quantum 
well system. This was not entirely expected as even a contact 
interaction still has contributions beyond on-site interaction 
especially for relatively small but finite barrier widths. In the 
latter case the Hubbard model reproduces at least the quali- 
tative trend, with a maximum discrepancy of ^ 20% for the 
parameters considered. We have compared the results from 
the Hubbard model also with the limiting case when no bar- 
rier exists between sites and a single QW is arbitrarily divided 
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into M equal sectors, each sector corresponding to one 'site'. 
Surprisingly, when enough 'wells' are considered, the Hub- 
bard model reproduces the entanglement within few percent. 
We interpreted this as the Hubbard model sites being a fine 
enough mesh discretization of the continuous spatial variable. 

We conjectured a theoretical maximum value for the aver- 
age single-site entanglement of two-particle trapped within M 
wells. We saw that the maximum value was not reached ex- 
cept in the case of 2 wells. We argued that for M > 2 some 
spatially modulated particle-particle interaction is needed to 
reach the maximum average single- site entanglement as to 
counteract the propensity of the particles to occupy the inner 
wells. 

Despite the calculated values of Uw and tw appearing to 
give very good results for relatively wide interwell barriers, 
we found that an even better match between the Hubbard 
model and the electron system could be achieved by rescaling 
the value of 11^/ {Cutw). This suggests that there were some 
small contributions to the interaction beyond the on-site repul- 
sion for the chosen well parameters, but that the main approxi- 



mation used — hopping parameters and interaction strength in- 
dependent of the site and estimated from the ground state of 
a single finite quantum well — ^remains valid. However, as the 
interwell barrier width decreases, these approximations fails 
and no rescaUng of Uw/ {Cutw) could improve the match be- 
tween the quantum well system and the Hubbard model re- 
sults. 

We also considered an attractive interaction [/ < and rel- 
atively wide interwell barriers. In this case the average single- 
site entanglement of the quantum well system was well ap- 
proximated by the Hubbard model. 

Finally we have considered a different type of 
entanglement — the spatial entanglement — for the quantum 
well system. Our results showed that the spatial entanglement 
tends to display in most parameter regions an opposite trend 
in respect to the average single-site entanglement. 

Future work includes considering long range Coulomb in- 
teractions, and how this affects comparison with the results 
from the Hubbard model. 
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